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We  find  an  algebraic  formula  for  the  A-partite  concurrence  of  N  qubits  in  an  X-matrix.  X-matrices 
are  density  matrices  whose  only  non-zero  elements  are  diagonal  or  anti-diagonal  when  written  in 
an  orthonormal  basis.  We  use  our  formula  to  study  the  dynamics  of  the  IV-partite  entanglement  of 
N  remote  qubits  in  generalized  A-party  Greenberger-Horne-Zeilinger  (GHZ)  states.  We  study  the 
case  when  each  qubit  interacts  with  a  local  amplitude  damping  channel.  It  is  shown  that  only  one 
type  of  GHZ  state  loses  its  entanglement  in  finite  time;  for  the  rest,  A-partite  entanglement  dies 
out  asymptotically.  Algebraic  formulas  for  the  entanglement  dynamics  are  given  in  both  cases.  We 
directly  confirm  that  the  half-life  of  the  entanglement  is  proportional  to  the  inverse  of  N .  When 
entanglement  vanishes  in  finite  time,  the  time  at  which  entanglement  vanishes  can  decrease  or 
increase  with  N  depending  on  the  initial  state.  In  the  macroscopic  limit,  this  time  is  independent 
of  the  initial  entanglement. 
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Even  though  entanglement  was  already  promoted  by 
Schrodinger  [1]  as  a  fundamental  aspect  of  quantum  the¬ 
ory,  while  in  mathematics  it  predated  quantum  mechan¬ 
ics  by  decades  [2] ,  its  value  as  a  resource  for  a  wide  range 
of  potential  applications  was  not  appreciated  until  re¬ 
cently  [3,  4].  Furthermore,  although  its  importance  is 
now  largely  recognized,  witnessing  and  quantifying  the 
entanglement  of  arbitrary  mixed  states  are  still  open 
questions.  In  only  few-party  cases  do  prescriptions  ex¬ 
ist  for  determining  the  entanglement  of  a  mixed  state 

[5-7]. 

The  problem  becomes  much  more  difficult  for  gen¬ 
uinely  multipartite  entanglement,  entanglement  shared 
between  more  than  two  parties.  Multipartite  entangle¬ 
ment  ( N  >  3)  is  thought  to  play  an  essential  role  in 
many  phenomena  including  quantum  metrology  [8]  and 
quantum  phase  transitions  [9].  Furthermore,  it  is  of  fun¬ 
damental  importance  to  understand  the  dynamics  of  mul¬ 
tipartite  entanglement  when  the  number  of  parties  shar¬ 
ing  entanglement  approaches  the  macroscopic  limit,  i.e. , 
N  — >  oo.  Previous  studies  of  the  dynamics  of  multipar¬ 
tite  entanglement  have  utilized  measures  that  fail  to  cap¬ 
ture  exactly  when  multipartite  entanglement  disappears 
[10-13]. 

An  essential  step  was  taken  by  Aolita  et  al.  [13]  where 
the  authors  utilized  the  entanglement  of  different  bipar¬ 
titions  of  an  X-qubit  system  to  qualitatively  study  the 
scaling  laws  for  the  decay  of  multiqubit  entanglement. 
The  caveat  to  this  approach  lies  in  the  fact  that  the  en¬ 
tanglement  of  different  bipartitions  is  a  necessary  but 
not  sufficient  condition  for  A-partite  entanglement.  The 
genuinely  multipartite  entanglement  between  X  parties 
can  vanish  before  the  entanglement  of  any  of  the  bipar¬ 
titions  vanish.  The  lack  of  such  analysis  is  mostly  due  to 
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the  fact  that  although  there  have  been  many  attempts 
to  solve  the  problem  of  determining  the  multipartite  en¬ 
tanglement  of  a  given  state  [14-16]  (see  also  references 
in  the  paper  by  Ma  et  al.  [17]),  an  algebraic  and/or  nu¬ 
merically  efficient  prescription  has  not  yet  emerged.  An 
algebraic  prescription  would  be  especially  desirable  since 
it  can  potentially  open  the  door  for  a  wide  range  of  ana¬ 
lytical  investigations  of  entanglement  dynamics. 

Recently,  based  on  previous  works  by  Pope  and  Mil- 
burn  [18]  and  Love  et  al.  [19],  a  new  measure  of  multi¬ 
partite  entanglement,  called  genuinely  multipartite  (GM) 
concurrence,  has  been  proposed,  and  it  has  been  shown 
that  GM  concurrence  is  an  entanglement  monotone  [17]. 
This  measure  reduces  to  Wootters’s  original  concurrence 
[6]  for  two  qubits.  Additionally,  an  algebraic  formula  for 
a  lower  bound  of  the  GM  concurrence  has  been  found 
by  Ma  et  al.  [17].  The  lower  bound,  when  calculated  for 
a  two-qubit  X-matrix,  matches  the  value  of  Wootters’s 
concurrence. 

The  X-matrix  of  Yu  and  Eberly  [20]  is  a  density  ma¬ 
trix  of  N  qubits,  written  in  an  orthonormal  product  ba¬ 
sis,  whose  non-zero  elements  are  only  diagonal  or  anti¬ 
diagonal.  The  concurrence  of  a  two-qubit  X-matrix  takes 
a  very  simple  form  [20]  and  that  is  why  these  two-qubit 
states  have  been  extensively  used  in  studying  the  dynam¬ 
ics  of  entanglement  between  two  qubits  in  many  scenarios 
[20-23], 

In  view  of  the  fact  that  the  GM  concurrence  lower 
bound,  derived  in  Ref.  [17],  matches  the  exact  value  of 
concurrence  for  a  two-qubit  X-matrix,  one  might  wonder 
if  the  lower  bound  might  also  be  exact  for  more  than  a 
two-qubit  X-matrix.  In  this  paper  we  will  prove  that  this 
conjecture  is  correct.  The  lower  bound  provided  by  Ma 
et  al.  [17]  is  realized  by  X-matrices.  We  thus  present  an 
algebraic  formula  for  the  GM  concurrence  of  an  X-qubit 
X-matrix.  This  is  our  principal  result  and  it  enables  ana¬ 
lytic  formulation  of  dynamics  of  X-partite  entanglement 
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in  different  scenarios.  For  illustration,  we  utilize  our  for¬ 
mula  to  directly  examine  the  decay  of  TV-qubit  entangle¬ 
ment  exposed  to  local  decoherence  channels. 

TV-partite  entanglement  is  defined  by  its  opposite, 
biseparability.  A  pure  TV-partite  system  \tjj)  S  Hi  ® 
H2  •  •  ■  ®  Hn  is  biseparable  if  there  is  a  bipartition  of  the 
TV  parties  Hi  G>  %2  •  ■  ■  ®Hn  =  Ha  ®  Hb ,  where  Ha  = 
Hj1®Hj2<8>-  ■  -®Hjk,  Hb  =  Hjk+1  <S)Hjk+2<S>-  ■  -®HjN  for 
which  | -ip}  =  \1>a)®\iI>b),  and  | ^a)  G  Ha  and  \i/jb)  G  Hb 
[3,  24].  In  other  words,  a  pure  state  is  biseparable  if  it  has 
at  least  one  pure  marginal  (reduced  density  matrix).  An 
TV-partite  state  that  cannot  be  written  as  an  ensemble  of 
biseparable  states  is  an  TV-partite  entangled  state. 

Before  introducing  GM  concurrence,  let  us  introduce 
the  set  of  all  bipartitions  of  TV  parties.  Each  biparti¬ 
tion  is  a  division  of  the  set  {1,2,...,  TV}  into  two  non¬ 
overlapping  and  non-empty  subsets.  The  set  of  all  such 
bipartitions  is  denoted  by  J  =  { Ji,  J2, . . .  J^n-i-i}. 
For  example  for  TV  =  3,  there  are  three  bipartitions 
Ji  =  {1|2,  3},  J2  =  {2|1, 3},  and  J3  =  {3|2,1},  so  that 
J  =  {Ji,  J2,  J; 3}- 

For  a  pure  state  | ip),  to  each  of  the  elements  of  J  we 
can  associate  two  reduced  density  matrices  A(\tp),Jj), 
and  B (| if)1  Jj)  by  tracing  out  either  of  the  subsystems 
associated  with  that  bipartition.  The  biseparability  of  a 
pure  state  can  be  determined  by  whether  for  any  of  the 
elements  of  J,  Jj)  and  B(\ip),  Jj)  are  pure.  If  so, 

|  i/j)  is  biseparable.  Thus  the  purity  of  the  jth  bipartition, 
to  be  denoted  by  11^(1^)),  is  a  key  parameter. 

For  a  pure  state  genuinely  multipartite  concurrence 
(GM  concurrence)  is  then  defined  [17]  as, 


Now  we  are  ready  to  extend  the  definition  of  GM  con¬ 
currence  to  all  mixed  states: 


Cgm(p)  =  min  Ca{p).  (3) 

a. 

If  Cgm  (p)  =  0  this  means  that  there  is  an  a  for  which 
Ca(p)  =  0.  Then  p  can  be  written  as  a  sum  of  pure  bisep¬ 
arable  states,  so  p  is  biseparable.  If  Cgm(p)  >  0,  there 
is  no  a  for  which  the  |V’f  )’s  are  all  biseparable  and  thus 
p  is  an  TV-partite  entangled  state.  The  GM  concurrence 
of  TV  parties,  as  defined  in  Ref.  [17],  is  a  monotone  of 
genuinely  multipartite  entanglement;  it  distinguishes  be¬ 
tween  biseparable  and  TV-partite  entangled  states,  is  con¬ 
vex,  invariant  under  local  unitary  transformations,  and 
non-increasing  under  local  operations  and  classical  com¬ 
munications  (LOCC)  [17].  The  operational  meaning  of 
GM  concurrence  in  terms  of  mutual  information  is  ex¬ 
plicitly  discussed  in  [26]. 

If  the  orthonormal  basis  for  the  X-matrix  is 
{|0, 0, . . . ,  0),  |0, 0, . . . ,  1), . . . ,  j  1, 1, . . . ,  1)},  then  one  can 
always  write  an  X-matrix  in  the  form  given  below 


/  ai 


X  = 


Ci2 


Zl  \ 


Z2 


biJ 


(4) 


Cgm{W)  ■  =  min  V2J1  -  B-j(\ip))- 
3  v 

Clearly,  Cgm(\4’))  >  0  and  it  is  equal  to  zero  if  and  only 
if  \t/j)  is  a  biseparable  state.  For  a  bipartite  system  this 
definition  reduces  to  the  I-concurrence  [25] . 

To  determine  whether  a  mixed  state  p  is  biseparable 
or  not  one  has  to  determine  whether  p  can  be  written 
as  a  convex  sum  of  pure  biseparable  states.  Thus  one 
has  to  check  all  the  ways  p  can  be  written  as  a  convex 
sum  of  pure  states  (all  pure  state  decompositions).  Let 
us  distinguish  different  pure  state  decompositions  of  p  by 
assigning  a  continuous  superscript,  a,  to  label  them: 

/3=E^“}  iv’fxci-  (i) 

i 

To  determine  whether  a  particular  pure  state  decompo¬ 
sition  is  a  sum  of  biseparable  states  or  not,  we  can  cal¬ 
culate  the  average  pure  state  GM  concurrence  for  that 
particular  a 

ca(p)  =  j2p(?)cGM(  m) 

i 

=  E  p(>a)  { m/n  ■  (2) 

7.  ^  ' 


where  n  =  2N  1,  and  we  require  \zi\  <  y/aibi  and  + 
b{)  =  1  to  ensure  that  X  is  positive  and  normalized.  One 
can  see  why  density  matrices  in  this  class  are  called  X- 
matrices.  It  can  be  shown  that  the  GM  concurrence  of 
an  TV-qubit  X-matrix  is  given  by 

Cgm  =  2  max{0,  \zi\  -  Wj},  i  =  0, 1, . . . ,  n  (5) 

where  Wi  =  Xq/j  V %  bj .  A  detailed  proof  of  this  result 
is  given  in  the  Appendix. 

Robustness  of  N-partite  entanglement:  Restricted 
forms  of  X-matrices  of  more  than  two  qubits  have  been 
used  in  some  recent  studies  of  the  dynamics  of  multi¬ 
partite  entanglement.  The  entanglement  measures  uti¬ 
lized  in  these  studies  yield  qualitative  information  about 
the  multipartite  entanglement  [27,  28].  Direct  study  of 
the  dynamics  of  genuinely  multipartite  entanglement  has 
been  an  elusive  problem  mainly  due  to  the  lack  of  an  an¬ 
alytical  measure  of  genuinely  multipartite  entanglement 
that  is  simple  to  calculate.  Only  for  Greenberger-Horne- 
Zcilinger  (GHZ)  states  that  undergo  pure  dephasing,  has 
there  been  a  successful  attempt  that  uses  a  geometric 
measure  to  give  the  exact  dynamics  of  TV-partite  entan¬ 
glement  [29].  Our  GM  concurrence  formula  provides  an 
opening  to  quantitatively  examine  the  conjectures  of  such 
studies  and  many  other  scenarios  whenever  the  initial 
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density  matrix  is  an  X-matrix  and  the  X  nature  of  the 
density  matrix  is  robust  in  the  dynamics. 

In  the  following  we  use  our  formula  to  directly  study 
the  dynamics  of  the  multiqubit  entanglement  shared  by 
N  qubits  when  each  qubit  is  subjected  to  a  local  am¬ 
plitude  damping  channel.  This  can  represent,  e.g.,  the 
spontaneous  decay  of  N  two-level  atoms,  each  in  a  sep¬ 
arate  zero-temperature  Markovian  reservoir.  For  a  two- 
level  atom  in  zero-temperature  Markovian  reservoir,  the 
evolution  of  ground  and  excited  states,  \g,  0)  and  |e,0), 
is  given  by 

U{t)\e,  0)  =  At\e,  0)  +  Bt\g,  1), 

U{t)\g,0)  =  |  <?,  0) ,  (6) 


FIG.  1.  The  versus  P  for  different  numbers  of  qubits  for 
the  initial  state  -|).  N  =  2  is  the  dotted  line,  N  =  10  is 
the  dashed  line  and  N  =  100  is  represented  by  the  solid  line. 


here  U{t)  is  the  local  propagator,  At  =  \/l  —  Pt ,  and 
Bt  =  VPf  Although  Pt  has  a  time  dependence  Pt  = 
1  —  e_7t  where  7  is  the  damping  rate,  we  can  also  think 
of  Pt  as  the  probability  of  the  decay,  and  write  everything 
as  a  function  of  P  instead  of  an  explicit  dependence  on 
time.  Thus,  in  the  following  we  drop  the  explicit  time 
dependence  of  P.  The  state  |1)  denotes  an  excited  state 
of  a  local  reservoir. 

We  study  the  dynamics  of  the  multiqubit  entanglement 
that  is  shared  initially  by  IV-partite  GHZ  states. 

\&P,a)  =cosa\e®N~kg®k)  +  sma\g®N~ke®k) 


In  Fig.  1,  we  plot  (/'ll1  versus  P  for  N  =  {2,10,100} 
qubits.  It  confirms  that  bulk  of  the  initial  entanglement 
dies  out  faster  (at  smaller  P’s)  as  the  number  of  qubits 
increases.  In  the  formula,  the  non-negative  factor, 
2(cos2a)  (1  —  P)^",  determines  the  decay  of  entangle¬ 
ment  for  N  2,  and  one  can  show  that  for  the  ampli¬ 
tude  damping  channel  the  half-life  of  the  entanglement 
depends  on  N  as 

Phalf-life  »  (9) 


This  is  a  GHZ  state  where  either  (IV  —  fc)’s  of  the  qubits 
are  initially  excited  and  the  rest  are  in  their  ground  state 
or  k  qubits  are  initially  excited  and  the  rest  in  their 
ground  states.  We  first  study  the  k  =  0  case.  We  present 
a  detailed  analysis  only  for  the  three  qubit  case  but  the 
generalization  to  N  qubits  is  straightforward.  By  tracing 
out  the  reservoirs  we  find  the  density  matrix  of  the  three 
atoms, 


/  ai 


/V°>  (*) = 


z  1 


zi  \ 


a  2 


02 


02 


bij 


(7) 


where 


ai  =  cos2  a|At|6, 
a,2  =  cos2  a\A^Bt\2. 
Z\  =  sin  a  cos  a  A ^ . 


b1  =  sin2  a  +  cos2  a|Pt|6, 
62  =  cos2  a\  AtB2\2, 


For  an  initial  |  ,  a)  state  the  concurrence  reads 

C^]  =  maxjOjQ^},  (8) 

=  2(cos2  a)  (1  —  P)"5"  tana|  —  (2iv_1  —  1)P^  J  . 


which  is  the  same  as  the  half-life  of  the  coherence  ele¬ 
ments  in  the  density  matrix.  We  observe  from  Fig.  1 
that  the  half-life  of  the  entanglement  decreases  as  the 
number  of  the  qubits  increases.  One  might  expect  a  sim¬ 
ilar  dependence  on  N  for  the  time  at  which  the  entangle¬ 
ment  disappears  completely.  On  the  contrary  this  is  not 
always  the  case.  In  order  to  show  this  we  solve  the  equa¬ 
tion  =  0  for  the  critical  value  of  P,  beyond  which 
the  concurrence  is  zero: 


P  > 


ltanQl  =  P 

2n~  i-P 


(10) 


If  Pc  <  1,  then  the  entanglement  has  a  finite  life  [22]. 
Otherwise  the  entanglement  dies  out  asymptotically.  In 
Fig.  2,  we  plot  Pc  versus  the  number  of  the  qubits  for  dif¬ 
ferent  initial  states.  We  observe  that  the  critical  value, 
Pc,  can  increase,  decrease  or  even  decrease  and  then  in¬ 
crease  as  a  function  of  N.  The  parameter  that  determines 
this  peculiar  dependence  of  Pc  on  N  is  tan  a,  which  one 
can  think  of  as  a  distance  of  the  initial  state  to  the  final 
state.  In  the  macroscopic  limit,  N  — >  00,  even  this  depen¬ 
dence  on  tan  a  is  suppressed.  Thus  although  the  half-life 
of  macroscopic  entanglement  is  very  small,  a  non-zero 
entanglement  lasts  for  a  constant  interval  of  time  before 
vanishing  completely. 

Similar  unusual  behaviors  were  observed  by  Aolita 
et  al.  [13]  for  the  entanglement  of  different  bipartitions  of 
the  N  qubits.  They  had  derived  similar  IV-dependence 
for  the  half-life  and  also  provided  examples  of  initial 
states  giving  Pc  increasing  with  N.  It  should  be  pointed 
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out  that  since  the  entanglement  of  different  bipartitions 
is  not  a  sufficient  condition  for  N-partite  entanglement 
similar  behavior  was  not  a  foregone  conclusion. 

For  the  k  >  0  case  one  can  show  that  the  initial  entan¬ 
glement  only  decays  asymptotically.  Below  we  present 
the  argument  for  TV  =  3  and  k  =  1  but  the  extension 
to  higher  TV  is  straightforward.  The  asymptotic  decay  of 
entanglement  is  due  to  the  fact  that  (e,  e,  e|p$(i)  |e,  e,  e), 
(e,ff,e|p$(i)|e,5,e),  and  (g,e,e\p^i)\g,e,e)  remain  zero 
for  all  times.  To  show  this,  we  note  that  our  initial  state 
is  a  superposition  of  two  possibilities.  Either  two  of  the 
atoms,  {1,2},  are  excited  and  the  other  atom,  {3},  is 
in  vacuum  state,  or  that  single  atom  is  excited  and  the 
atoms  {1,2}  are  in  their  ground  states.  Since  all  reser¬ 
voirs  are  initially  in  their  ground  states,  if  an  atom  is  in 
its  ground  state  initially  it  will  always  remain  there.  But 
the  three  diagonal  terms  that  we  referred  to  require  the 
atom  {3}  and  at  least  one  of  the  other  two  atoms  to  be  ex¬ 
cited  simultaneously.  Since  this  is  forbidden,  all  of  these 
matrix  elements  remain  zero.  Thus  the  negative  con¬ 
tribution  to  the  concurrence  formula  remain  zero.  This 
argument  can  be  generalized  for  TV  >  4  to  all  of  the  GHZ 
states  except  for  |$^,a),  because  in  the  state 

all  of  the  atoms  can  be  initially  excited.  The  concurrence 
of  all  |  ,  a)  initial  states  with  k  >  0  is  given  by 

C$  =  |sin2a|(l-P)£. 

Thus  for  k  >  0,  concurrence  only  dies  when  P  =  1. 

N-partite  entanglement,  either  as  a  resource  for  quan¬ 
tum  computation  or  as  a  fundamental  property  of  quan¬ 
tum  theory,  has  been  difficult  to  quantify,  especially 
for  mixed  states.  This  is  an  important  drawback  since 
many  of  the  algorithms  in  quantum  computation  need  N- 
partite  entanglement  between  a  large  number  of  qubits, 
and  inevitable  interaction  of  these  qubits  with  the  envi¬ 
ronment  renders  initial  pure  states  mixed  and  diminishes 
their  entanglement.  Thus  it  is  of  interest  to  understand, 
quantitatively,  the  dynamics  of  N-partite  entanglement 
when  the  qubits  sharing  it  come  in  contact  with  different 
environments.  Here,  we  have  found  an  algebraic  formula 
for  the  genuinely  multipartite  (GM)  concurrence  of  N- 


FIG.  2.  Pc  versus  number  of  qubits  for  different  initial  states 
|<F^,a).  From  bottom  up  tana  =  0.01,0.1,0.2,0.5  and  1 
respectively. 


qubit  density  matrices  that  can  be  written  as  X-matrices 
in  an  orthonormal  product  basis.  This  development  al¬ 
lows  TV-partite  entanglement  to  be  quantified  for  such 
states.  The  formula  opens  up  the  possibility  of  study¬ 
ing  entanglement  dynamics  of  TV-qubit  states  in  different 
scenarios,  as  long  as  the  X-form  of  the  density  matrix  is 
preserved. 

Using  the  concurrence  formula,  we  have  studied  the  dy¬ 
namics  of  TV-partite  entanglement  of  TV  two-level  atoms 
interacting  with  local  amplitude  damping  channels.  We 
showed  that  only  for  k  =  0  the  |<b^,a)  initial  states 
lose  their  N-partite  entanglement  in  finite  time.  Alge¬ 
braic  formulas  for  the  concurrence  were  presented.  It  is 
observed  that  for  large  TV  the  bulk  of  initial  concurrence 
decays  with  a  rate  inverse  to  TV.  For  a  given  TV  and 
k  =  0,  the  time  at  which  entanglement  vanishes  to  zero 
is  determined  by  the  distance  of  the  initial  state  from  the 
final  state.  In  the  macroscopic  limit  this  time  is  indepen¬ 
dent  of  a  too.  An  open  question  is  whether  this  time 
interval  also  appears  for  other  kinds  of  initial  states  in 
the  macroscopic  limit  and  whether  it  has  any  observable 
effect. 
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I.  APPENDIX:  CONCURRENCE  OF  TV-QUBIT 
X-STATES 

In  [17],  Ma  et  al.  presented  a  lower  bound  for  the  GM 
concurrence.  The  lower  bound  of  GM  concurrence,  de¬ 
rived  in  [17],  for  an  X-matrix  is  given  by 

Cgm  >  2  max{0,  \zf\  -  wt},  *  =  0,  l,...,n  (11) 

where  Wi  =  YZ'j^i  \J aYair  In  thc  following  we  will  show 
that  this  lower  bound  is  exact  for  all  X-matrices.  With¬ 
out  loss  of  generality  we  can  assume  that  y/a\b\  >  y/aibi- 
Since  we  have  assumed  that  y/a\b\  >  y/aJh,  it  is  easy  to 
show  that  \zi\  —  Wi  <  0  for  i  >  1,  so  that  Eq.  11  reduces 
to 

Cgm(X)  >  2  max{0,  \z\\  —  u>i}  (12) 

We  will  show  that  this  bound  is  actually  an  equality. 
First  let  us  prove  a  lemma  that  we  will  utilize  in  our 
proof. 

Lemma  1.  The  GM  concurrence  of  an  X-matrix  for 
which  a\bi  >  cq&i,  and  aj  =  bj  =  0  for  all  j  ^  {*,  1}, 
is 

Cgm  (Xu)  =  2  max{0,  |zi|  -  \/aA}  (13) 

Proof.  We  already  know  that  this  quantity  is  a  lower 
bound  of  GM  concurrence.  Thus  we  only  need  to  show 
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that  it  is  also  an  upper  bound.  We  will  do  this  by  map¬ 
ping  Xu  to  a  two-qubit  density  matrix,  R,  and  then  show 
that  CGM(XU)  is  bounded  from  above  by  Wootters’s 
concurrence  of  R,  where  C(R)  =  2max{0,  \zi\  —  Vaib,}. 

Before  going  forward  let  us  introduce  some  notation. 
Since  we  are  working  with  qubits,  we  can  represent  each 
vector  (ket)  of  the  above  basis  as  a  number  from  0 
to  2N  —  1  written  in  the  binary  basis.  For  example, 
|0, 0,  •  •  •  ,  0)  =  |0),  |0,  -  -  -  ,0,1)  =  |1),  |0,  ■  ■  -  ,1,0)  =  |2), 
and  so  on  |1,  ■  ■  •  ,  1, 1)  =  \2N  —  1).  We  also  denote  the 
bit-flipped  states  in  the  same  way,  |i)  =  \2N  —  i  —  1), 
e.g.,  |0)  =  12^  —  1).  In  places  where  we  need  to  label  the 
individual  qubits,  we  will  do  so  by  using  a  subscript  on 
the  bits, 

We  perform  the  mapping  by  focusing  on  a  specific  bi¬ 
partition  of  the  qubits.  The  four  non-zero  diagonal  ele¬ 
ments  of  Xu  are  {ai,  a*,  bi,  bi },  corresponding  to  projec¬ 
tors  {|0)(0|,  \i—  l)(i  — 1|,  1 0) (0 1 ,  and  \i  —  1  )(i—  1 1}  respec¬ 
tively.  Those  qubits  that  contribute  1  to  the  ket  \i—  1)  we 
designate  as  party  F.  The  rest  of  the  qubits  we  denote 
as  party  G.  For  example,  with  7  qubits,  which  we  denote 
as  (91, 92, 93, 94, 95, 96, 97),  where  i  =  6,  the  basis  states 
are 

|0)  =|0i,02,03,  04,0b,  08, 07), 

1 5)  =|0i,02,03, 04, 15, 06, 17), 

1 127)  =  1 1 1 , 1 2 , 13, 14, 15, 16, 17), 

1 122)  =  1 1 1 , 12, 13, 14,  05, 16, 07).  (14) 


two-qubit  system  back  to  the  multi-qubit  basis  states. 
We  pick  the  PSD  whose  average  concurrence  is  the  min¬ 
imum  amongst  all  possible  PSD’s  of  R.  Thus 

R  =  Y^Pi\i>i){^i\,  C(R)='£Pic(Wi))-  (17) 

i  i 

In  Eq.  17,  C(R)  is  Wootters’s  concurrence,  which,  by 
definition,  is  equal  to  the  minimum  average  concurrence 
over  all  possible  PSD’s  of  R.  As  mentioned  before,  each 
pure  state  | ipi)  can  be  mapped  back  to  an  N-qubit  state 
(h pi)  |4>i)),  producing  a  PSD  for  Xu, 

Xli=Y,Pi\*i)(*i\-  (18) 

i 

Since  GM  concurrence  is  convex  by  definition,  we  have 

Com <  J^P^GM^i))-  (19) 

i 

For  a  pure  state  the  GM  concurrence  is  defined  by 

CGM(\*i))  =  min^/l-n.d^)),  (20) 

3  v 

where  the  minimum  is  taken  over  all  bipartitions,  J,  of 
the  N  qubits.  Therefore,  the  GM  concurrence  must  be 
bounded  by  any  specific  bipartition,  including  the  bipar¬ 
tition  of  the  N  qubits  to  party  F  and  party  G. 


Then  party  F  is  given  by  qubits  (91, 92, 93,  94, 96),  and  the 
remaining  two  qubits,  (95, 97),  make  party  G.  Under  this 
bipartition,  we  can  write  Xu  using  the  following  basis 
states, 

I  If)  =  |0i,  02,  03, 04, 06),  |  1>)  =  |li,  I2, 13, 14, 16), 

I  4-g)  =  {IO5,  O7),  |  Tg)  =  II5, 17), 

Xu  =  ai|  4-f4-g)  (4-fJ-g  I  +  bi\  1>tG)(1>tG  I 

+  Oil  4-FtG)(4-FtG  I  +  bi\  t F-I-g) (t f4-G  ] 

+  Zi\  -J-F  4-g) (t Ft G  |  +  Zl\  tFtG)(4-F4-G  I 

+  Zi\  |FtG)(tFlG  |  +  Z*  |  IfIgX-I-fIg  |,  (15) 


Cgm( |*4))  <  v^i-n^d^)).  (21) 

Using  the  same  mapping  as  done  for  Eq.  16  it  is  easy 
to  show  that  C'(|V’i))  is  equal  to  the  right  hand  side  of 
Eq.  21.  Therefore  we  conclude  that 

Cgm(Xu)  <  Y,PiCGM(\*i))  <  J2P'C (1^)) 

i  i 

=  C(R)  =  2max{0,  \zi\  -  y/dibi},  (22) 

where  the  right  most  equality  is  found  by  evaluat¬ 
ing  Wootters’s  concurrence  for  R  under  the  assumption 
ai&i  >  Oi&j.  This  upper  bound  matches  the  lower  bound 
and  therefore  it  is  the  exact  value  of  Cgm(Xu)-  □ 


We  see  that  if  we  restrict  attention  to  the  subspace  de¬ 
fined  by  the  non-zero  elements  of  Xu  we  can  map  Xu 
to  a  two  qubit  density  matrix,  R,  which,  in  the  basis 
{|  IfIg),  I  J-fIg),  |  1>1g),  |  tFtG)},  reads 


Xu 


R  = 


CLi 

* 

Z 


(16) 


Now  that  we  have  a  two-qubit  density  matrix,  we  can 
take  advantage  of  Wootters’s  concurrence.  Note  that 
from  each  pure-state  decomposition  (PSD)  of  R  one  can 
make  a  PSD  of  Xu  by  mapping  the  basis  states  of  the 


Next,  we  generalize  this  result  to  all  the  X-matrices. 
We  do  so  by  decomposing  the  X-matrix  into  a  convex 
sum  of  Xu  matrices.  Let  us  first  look  at  the  case  for 
which  \zi  |  —  wi  >  0. 

(a)  \zi\  —  W\  >  0.  Note  that  \zi\  —  Wi  <  0  for  i  >  2 
since  y/a±bi  >  y/o-ibi  >  \zi\.  First  by  a  change  of  phase 
of  the  basis,  which  is  a  local  unitary  transformation,  we 
change  z\  to  \zi\.  This  only  changes  the  phase  of  the 
other  off-diagonal  elements.  Then  we  decompose  X  in 
the  following  form. 


X  =  A  +  £  Si, 

i>  1 


(23) 
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where  Si  is  an  Xu  matrix  whose  two-qubit  counterpart 
reads 


Ri  = 


^  %%  V 

Cli  Zi 
Z*  bi 

\  Vaibi  Vi  ) 


(24) 


where 


.  ^  w\  .  .  .  w i  . 

^11  —  CL  i(l - ),  A2n,2n  =  0  l(l -  ), 

V  CL\b\  Y  CL\b\ 

Alt2  n  =  ^2n,l  =  \Zl\  —  tUl, 

=  0  *^{l,2n},or  j^{l,2n}, 


Xi  = 


CL  1 V  CLiOi 
W  i 


Z/i  = 


&i\/a 


1 

u>i 


(25) 


It  can  be  shown  that  S)’s  are  all  proportional  to  valid 
density  matrices,  since  they  are  non-negative  hermitian 
matrices.  The  proportionality  constant  is  between  zero 
and  one  and  can  be  interpreted  as  probability.  Us¬ 
ing  Lemma  1  one  can  show  that  all  Si s  are  bisepara- 
ble  matrices  (though  not  normalized).  Regarding  the 
first  matrix  in  the  decomposition,  the  proportionality 
constant  is  An  +  A2Hi2n,  and  its  GM  concurrence  is 
2(|,zi|  —  wi)/{An  +  A2n,2n)-  Due  to  the  convexity  of 
GM  concurrence  we  conclude  that 

(An  +  A2n,2n)  |?(N  ~  Wl)  =  2(M  -  Wl)  (26) 

All  +  A2n,2n 


is  an  upper  bound  for  the  GM  concurrence  of  X.  Since 
2(|zi|  —w{)  is  also  a  lower  bound  for  the  concurrence,  it 
is  the  exact  value  of  the  GM  concurrence. 

Note  that  for  the  above  decomposition  to  work  we  had 
to  assume  that  \z\\  >  w\.  We  now  turn  to  the  case 
\zi\  <  wi .  We  seek  to  show  that  all  such  density  matrices 
are  biseparable.  We  consider  two  different  scenarios. 

(b)  y/aibi  >  wi.  In  this  case  the  matrix  X  can  be 
decomposed  to  matrices  similar  to  the  previous  case. 

n 

X  =  '52%,  (27) 

i>  1 


where  is  an  Xu  matrix  whose  two-qubit  counterpart 
reads 


/  ai  Ti 

ZiTi 

ai  Zi 

z*  h 

V  z\Ti 

biR 

(28) 


Since  |zi|T(  <  yj mbi  and  \zi\  <  y/aibiTi,  we  can  invoke 
Lemma  1  to  confirm  R'  is  biseparable  for  all  i.  The 
fact  that  R[  is  not  normalized  does  not  interfere  with 


the  proof  of  biseparability  as  one  can  always  factor  out 
Tr[R'].  Now  we  focus  on  the  last  case. 

(c)  y/aibi  <  wi .  In  this  case  we  divide  our  matrix 
into  two  positive  semi-definite  matrices  X  =  Ki(t,r)  + 
Ki(r,t). 


Zit  \ 

z2r 

nr 


b2r 

bit  / 


r  =  1  —  t. 

Note  that  since  wi  <  3y/aibi,  then  |  >  t  >  r.  One 
can  show  that 

t^/aibi  =  rwi,  (29) 

which  guarantees  that  A'i  (t,  r )  falls  in  the  category  of 
case  (a).  Since 

t\zi\  <  rw i, 

r\zj\  <  (t  —  r)y/aibi  +  rwj ,  (30) 

Ki(t,r)  is  biseparable.  Regarding  matrix  Ki(r,t),  since 
r\zi\  <  twi, 

n 

t\zj\  <  ry/aibi  +  t  Y  Aa.  (sc 

i¥-  i,i 


/  ait 


Ki(t,r)  = 


a2r 


anr  z 
z*r  b 


V  ztt 


where 


t  = 


w  1 


wi  +  y/a\bi 


and 


it  does  not  fall  in  the  category  of  case  (a)  and  thus  be¬ 
longs  to  either  case  (b)  or  case  (c).  If  it  falls  in  the 
category  of  case  (b)  then  we  can  conclude  that  it  is 
biseparable.  If  not,  we  divide  Ki(t,r)  into  two  matri¬ 
ces  Ki(t,r)  =  K2(t',  r')  +  K2(r',  t'),  as  before.  Each  time 
we  divide  a  matrix  in  this  way  the  trace  of  the  remaining 
part  is  strictly  smaller  than  the  trace  of  the  step  before: 
Tr [Ki(rl,tl)\  <  0.75L  Thus,  we  can  write  the  matrix  X 
as  a  convex  sum  of  biseparable  states  and  a  remaining 
part  that  can  be  made  arbitrarily  close  to  zero.  There¬ 
fore  matrix  X  is  a  biseparable  matrix.  This  completes 
the  proof  for  all  X-matrices.  Therefore,  we  have  proved 
that  the  GM  concurrence  of  a  iV-qubit  X-matrix  is 

Cgm{X)  =  2  max{0,  \zi\  -  i uj,  1  <  i  <  n  (32) 
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